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Abstract This paper investigates the potential for 
stabilizing an inverted pendulum without electric de¬ 
vices, using gravitational potential energy. We propose 
a wheeled mechanism on a slope, specifically, a wheeled 
double pendulum, whose second pendulum transforms 
gravity force into braking force that acts on the wheel. 
In this paper, we derive steady-state equations of this 
system and conduct nonlinear analysis to obtain pa¬ 
rameter conditions under which the standing position 
of the first pendulum becomes asymptotically stable. 
In this asymptotically stable condition, the proposed 
mechanism descends the slope in a stable standing po¬ 
sition, while dissipating gravitational potential energy 
via the brake mechanism. By numerically continuing 
the stability limits in the parameter space, we find that 
the stable parameter region is simply connected. This 
implies that the proposed mechanism can be robust 
against errors in parameter setting. 

Keywords Wheeled inverted pendulum • Passive 
control • Gravity • Friction • Asymptotic stabilization 


1 Introduction 

Electric and electronic control devices are indispens¬ 
able for a variety of modern technologies. However, 
these technologies typically become useless during mas¬ 
sive power outages such as those caused by natural 
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or other disasters. In this paper, we consider a non- 
electrified alternative control method to stabilize an 
inverted pendulum using gravitational potential. Our 
proposed mechanism is a wheeled inverted pendulum 
that descends a slope. The brake mechanism of our pro¬ 
posed mechanism transforms gravity force into friction 
force between the pendulum and the wheel. This friction 
produces a restoring force by which the pendulum is 
asymptotically stabilized in a standing position. 

Approaches similar to ours can be found in the field 
of passive dynamic walking, pioneered by McGeer [1], 
in which two-legged mechanisms are designed to sta¬ 
bly walk down a slope that consume only gravitational 
potential energy. Extensive studies have been reported 
on the passive dynamic walking, including experimen¬ 
tal development of passive walkers [1-3] and nonlinear 
analyses of passive dynamic walking based on simplified 
models [4-6]. Early insights into the use of such passivity 
can also be found in the study of passive gravity-gradient 
attitude stabilization [7-10] wherein the alignment of 
one axis of a satellite along the earth’s local vertical 
direction was achieved without the use of active control 
elements. 

On the other hand, the wheeled inverted pendulum 
has attracted significant attention in the fields of control 
engineering and robotics. Because of the applications 
of wheeled inverted pendulums in personal mobility de¬ 
vices, including the Segway® [11], methods for control¬ 
ling wheeled inverted pendulums have been developed 
via approaches such as partial feedback linearization 
[12], inclined surface control [13], sliding-mode velocity 
control [14], neuro-fuzzy-based control [15], and robust 
control based on a quasi-linear parameter-varying model 
[16]. Not surprisingly, these studies implied the use of 
electric devices. 
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Fig. 1 Friction-controlled wheeled inverted pendulum 
(FCWIP). 


In this paper, we merge concepts from passive dy¬ 
namic walking and studies of the wheeled inverted pen¬ 
dulum to derive our new mechanical design for the non- 
electrified stabilization of a wheeled inverted pendulum. 
As stated above, we propose a wheeled double pendulum 
mechanism, whose second pendulum transforms gravity 
force into braking force that acts on the wheel. To in¬ 
vestigate the dynamic stabilities of this newly proposed 
mechanism, we start with deriving a nonlinear analyt¬ 
ical model of the mechanism to examine the stabilities 
of its steady states. Three types of critical points arise 
in the analytical model. These critical points are ana¬ 
lytically characterized and numerically continued in the 
parameter space to obtain stability limits for the steady 
standing motions. It is found that the stability of the 
proposed mechanism is limited by Hopf bifurcation and 
vanishing external resistance on the wheel. 

2 Wheeled inverted pendulum with friction 
control 

We propose a wheeled inverted pendulum with friction 
control, as shown in Fig. 1, that comprises 1) a wheel 
placed on a slope without slipping, 2 ) a double pen¬ 
dulum suspended on the wheel axis, and 3) a friction 
control mechanism that generates a braking force on 
the wheel proportional to the angle between the first 
and second pendulums. Hereinafter, we refer to this 
model as a friction-controlled wheeled inverted pendu¬ 
lum (FCWIP). 

Configuration of the FCWIP can be described by 
a three-dimensional generalized coordinate: denotes 
transpose) 

e = {e^,02,e3f, (i) 


Table 1 Mechanical parameters of the wheeled double pen¬ 
dulum. 


Parameters 

Values 

mi 

mass of wheel 

0.1 kg 

m 2 

mass of 1st pendulum 

0.2 kg 

mz 

mass of 2nd pendulum 

1 kg 

h 

moment of inertia of wheel 

mir^/2 kg-m^ 

^2 

moment of inertia of 1st pendu¬ 
lum 

m2pfl2 kg-m^ 

h 

moment of inertia of 2nd pendu¬ 
lum 

0.25 kg-m^ 

r 

radius of wheel 

0.2 m 

1 

length of 1st pendulum 

1 m 

Ig 

placement of center of gravity of 
1st pendulum 

1/2 

w 

length of 2nd pendulum 

0.7 m 

WG 

placement of center of gravity of 
2nd pendulum 

0.6 m 

9 

acceleration of gravity 

9.8 m/s^ 

a 

angle of slope 

0.1 rad 


where 9i is the rotational angle of the wheel, 02 is the 
absolute slant angle of the first pendulum, and 03 is 
the relative angle of the second pendulum from the first 
pendulum. In addition, we consider the corresponding 
generalized force: 

T = (Ti,T2,r3)^, ( 2 ) 

where Ti is a torque acting on 0i. 

2.1 Wheeled double pendulum 

Unless the friction control mechanism (or T) is speci¬ 
fied, the FCWIP in Fig. 1 is simply a wheeled double 
pendulum whose Lagrangian is given by 

L:=T-U: (3) 

T := —0i(Qir + Ii) -\- —cos(03) -I- Q 2 + I 2 } 

+ + h) + MsQsiwG - I cos ( 6 » 3 )} 

+ 0 i^ 2 {-( 53 f cos(a - 02 - 03 ) + Qat cos(q; - 02 )} 
- cos(a - 6*2 - 6 * 3 ), (4) 

U := -g{Q 3 cos( 6 l 2 + 03 ) - Qi cos( 6 » 2 ) + 0iQi sin(a)}, 

( 5 ) 

with 

Qi = {mi + TO2 + TO3)r, 

Q 2 = m 2 lG + m3{wQ + 

Q3=m3WG, Q4 = m2lG + m3l, (6) 

where the physical parameters are listed in Table 1. 
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Fig. 2 Example of the structure of the brake mechanism 


On substituting L into Lagrange’s equations with 
the generalized force T, we obtain the equations of the 
motion of the wheeled double pendulum as 

M0 = F + T, (7) 

with 

= M, Mn = Qir + Ii, 

Mi 2 = -Qsr cos(a -62- O3) + Q^r cos(a - O2), 

Mi 3 = -Qsrcosia -62- ^3), 

M22 = COs(d3) + Q2 + 72) 

M23 = QsiwG — I cos(6i3)}, M33 = Q3WG + h, (8) 

and 

Fi =Q3r{02 + O3Y sin(a -62- O3) 

— QatO"^ sin(Q; — 6*2) + gQi sin(a), 

F 2 = — <33^3(202 + ^ 3 ) sin( 03 ) 

- g{Q 3 sin(6»2 + 6*3) - Q4 sin(6»2)}, 

F3 =Q3l0l sin(6»3) - gQ3 sin(6»2 + ^s), ( 9 ) 

where represents the (z, j ) component of the matrix 

M and Fi is the zth component of the vector F. 


2.2 Friction control mechanism 

Next, we introduce a friction control mechanism (FCM) 
into the wheeled double pendulum by specifying T as 
follows. 

Let 2 be a displacement of the brake rod outputted 
from the cam mechanism, as shown in Fig. 1, and sup¬ 
pose that the cam function 2 ( 03 ) is given as a linear 
function: 

2 = z (6<3) := p (6»3 - 77), z = p 93 (10) 

where p is a cam ratio and zy > 0 is an offset angle. 
Accordingly, the follower is expected to follow both the 
positive and negative rotation of the cam. 

Then, we consider a brake mechanism, as shown 
in Fig. 2. In this mechanism, a pad (in light gray) is 


bonded on a brake rod (in dark gray) and sandwiched 
between the fixed base on the first pendulum and the 
brake disk without clearance at z = 0. We refer to the 
lower half of the pad as a brake pad and the upper half 
as a dummy. The dummy pad has no function in terms 
of braking but is assumed to have the same mechanical 
property as the brake pad. 

Thus, the brake pad touches the brake disk when 
z > 0 but is separated from the disk when z < 0. We 
assume that the brake rod receives a continuous reaction 
force from the pads in the following form: 

R = R{z, z) := — {kbZ + Cbz) 

= - P {^ 6 ( 6*3 - zy)-b Ch^l, ( 11 ) 

where C{, and kb are viscoelastic coefficients of the pad. 

The reaction force R produces a torque on 63 as a 
generalized force Fg^ on 6 * 3 , given by 

Oz 

Ff)3 = ^)- (12) 

At the same time, we assume that R causes a Coulomb 
friction force between the brake pad and the brake disk. 
This force can be modeled by a tangential force on the 
contact surface as 

Fr = /zi?(z, z) sgn(0i - 02 )xiz) (13) 

where p, is the Coulomb friction coefhcient, sgn(-) is the 
unit signum function, and x(-) is the unit step function 
representing the separation of the brake pad from the 
brake disk. We have the torques R on 9i (z = 1, 2, 3) as 


(T3=rbFR-ci\93\9i, 

< T 2 = -uFr, (14) 

[r 3 = Fg^ = pi?(z,z), 

where ci is the coefficient of the quadratic resistance 
including aerodynamic force on the wheel (or 9i). Table 2 
summarizes the parameters of the FCM and quadratic 
resistance. Note that the value of the spring coefficient 
listed in Table 2 can be obtained approximately from a 
medium-carbon steel rod (Young’s modulus 205 GPa) 
of 5 X 10“"^ m diameter and 0.5 m length. 

Therefore, we derive the dynamic model of the 
FCWIP as the wheeled double pendulum in (7) with 
the braking torque in (14). 


2.3 Numerical examples 

Figure 3 shows a numerical solution of the FCWIP 
model obtained by solving (2), (7), and (14) from the 
trivial initial state Pi( 0 ) = 6 * 2 ( 0 ) = 0 , 6 * 3 ( 0 ) = rj (or 
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Table 2 Parameters of the FCM. 


Parameters 

Values 

rb 

radius of brake disk 

0.18 m 

P 

cam ratio 

1/20 

kb 

spring coefficient of the brake 

8 X 10-^ N/m 

Cb 

viscous coefficient of the brake 

(2 X lO"* Ns/m) 

P 

Coulomb friction coefficient of 
the brake 

(0.249) 

V 

offset of the brake mechanism 

(2 X 10“"* rad) 

Cl 

coefficient of quadratic resistance 
on the wheel 

(5 X 10—* Nms2) 


Parentheses around values denote nominal values. 



Fig. 3 Time responses of the FCWIP for the condition listed 
in Tables 1 and 2. 
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1,2,3). The area hatched in gray is the set of initial 
angles from which the state converges to fallen positions 
of the FCWIP model, and the white area surrounded 
by the hatched area is the set that converges to the 
steady standing state, as shown in Fig. 3. From Fig. 4, 
it appears that the set of the initial angles that belongs 
to the standing position forms a mostly connected area. 
Therefore, it can be expected that the proposed FCWIP 
model exhibits some robustness against disturbance of 
the initial angles. 

Considering potential real-life applications of the 
proposed system, dependencies of the parameters on 
the size of the attraction basin that belongs to the 
standing position could be a crucial problem. This will 
be addressed in a future study. 


3 Steady-state analysis 

3.1 State-space representation 

For a simple expression, we perform a time scale trans¬ 
formation t := qt*, where g is a time scale and t* is 
nondimensional time. Taking a state vector: 

X := {01,02,03 - V,Q&i,Q02,qO3)'^, 0i ■= d0i/dt, (15) 

we transform the dynamic model (7) and (14) to a 
state-space form: 

X = M{x)~^^f{x) + t{x)'^, x:=dx/dt*, (16) 

with 


Fig. 4 Attraction basin of the stable steady state in Fig. 3. 


M{x) := dia.g(^E 3 ,M{x)j, 


(17) 


2 ;( 0 ) = 0), and 0i(O) = 0 (i = 1,2,3). The parameter 
values are listed in Table 1 and 2 that are empirically cho¬ 
sen to achieve a stable standing motion. For numerical 
integration, a fourth-order Runge-Kutta-Gill method is 
employed with a time step of 10 “^. 

As shown in Fig. 3, the FCWIP model becomes 
asymptotically stable under the suitable parameter con¬ 
dition. In this example, the angle of the first pendulum 
02 (t) converges to a small negative value representing 
a standing position that is slightly slanted toward the 
upside of the slope. Consequently, the angle of the sec¬ 
ond pendulum 03 (t) converges to a small positive value 
that represents it hanging from the first pendulum to 
produce the brake force Fr in (13). Moreover, the de¬ 
scent velocity of the wheel 0i{t) converges to 9.46 rad/s 
(6.81 km/h in translational velocity). 

Figure 4 shows the sets of initial angles 02 ( 0 ), 03 ( 0 ) 
of the first and second pendulums, respectively. The 
other initial values are set as 0i(O) = 0i(O) = 0 {i = 


( Xi \ 


( ^ \ 

X 5 

, t(x) := 

0 

xe 

0 



\q^T{x)j 


where M{x), F{x), andT(a:) are the matrix and vectors 
in the dynamic model (7) and (14) via (15), and E3 is 
a 3 X 3 identity matrix. 

We choose q := (fc&p^)“^/^ to normalize the spring 
coefficient kb and introduce nondimensional parameters 
listed in Table 3. In this case, the components of the 
vectors in (18) are obtained as 

q'^Fi = -Q 3 r{x 3 + xqY sin (?7 + X2+X3- a) 

+QAf'x\ sin(a ;2 — a) -I- g*Qi sin(Q;), 

< q^F2 = -Qslxei^Xb + xq) sin (?7 -|- 0 : 3 ) 

-g*{Q3 sin(r? + X2 + X3) - Qi sin(x 2 )}, 

= Q3IXI sin (?7 -b X3) - g*Q 3 sin (?7 + X2 + X3), 
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Table 3 Nondimensional parameters. 


Parameters 

Values 

Definition 

9 

time scale 

200“^/^ 


K 

spring coefficient 

1 

{q^p^)kb 

9 * 

acceleration of gravity 

4.9 X 10“2 

(9^)9 

9 * 

Coulomb friction coeffi¬ 
cient 

(0.8964) 

{Vbp ^)/r 

F, 

viscous coefficient 

(3.536) 

iqp‘^)cb 


The values are transformed from the original parameters in 
Tables 1 and 2. Parentheses around values denote nominal 
values. 


and 

q^Ti = -fj,*sgn(x4 - X5)xix3) {xa + cIxq} 

-Ci\x4\x4, 

q^T2 = ^*sgn(x4 - X5)x(a;3) {xs + c^xq} , 
q'^Ts = - {xa + c^xe} . 


(19) 


( 20 ) 


3.2 Assumption of steady state 

In view of the numerical example presented in Fig. 3, 
we consider a steady state x* of the FCWIP model (16) 
that satisfies 

a; = jc* := (cc, 0 , 0 , 0 , 0 , 0 )^ , w > 0 (constant). ( 21 ) 

The steady state x* mentioned above describes the 
rotation of the wheel down the slope with a constant 
angular velocity x^ = uj > 0 while the first and second 
pendulum maintain certain steady angles X 2 and 
respectively, with x^ = Xq = 0. The angles X 2 and X 3 
are assumed to satisfy the following conditions: 

(A) X 2 < 0: the first pendulum reaches a standing posi¬ 
tion (slightly) slanted to the upside of the slope. 

(B) Xg > 0 (or 2 ;* >0): due to (A), the second pendulum 
hangs from the first pendulum (due to gravity and 
while maintaining X 2 < 0 ) to produce the brake force 
Fr in (13). 

These conditions are required to stabilize the first pen¬ 
dulum in a standing position. Because, they guarantee 
existence of the braking force Fr caused by a negative 
clearance between the brake pad and disk, Xg > 0 (or 
2 ;* > 0 ), and that is mechanically caused by X 2 < 0 . 
Otherwise, the braking force vanishes, and the FCWIP 
becomes nothing more than an uncontrolled wheeled 
double pendulum that can never be stabilized around 
the standing position. 

In addition, note that the FCWIP model, in absence 
of the floor model of the slope, can theoretically have 


another stable steady state, a static equilibrium where 
the second pendulum is hanging down at rest at X 2 = 
TT -b £2 and X 3 = £3 for small £ 2 , £3 > 0 . 


3.3 Steady-state equation 

The steady-state equation is given by 

X* =Mix*)-^{fix*)+Tix*)}, (22) 

where the derivative x* is the constant vector already 
defined in ( 21 ) and x* is an unknown vector representing 
the steady state. Multiplying both the sides by M{x*), 
we obtain 


f{x*) + T{x*) = M{x*)x* = diag(A 3 ,M(a:*))a:* = x* 

(23) 

where the last equality is due to the zero components 
of X* := (w, 0,0,0, 0,0). Therefore, the steady-state 
condition is obtained as follows: 


X4 = 

= u;, 

= x*e 

= 0, 

(24) 

0 = 


+ q^Ti = 

9 *Qi sin(a) — ciw^ — fi*X3, 

(25) 

0 = 

q^F 2 

+ q^T 2 = 

- 9 *{Q 3 sin (?7 + X2+ X3) 





- Q 4 sin(x2)} -b 9 *x'^, 

(26) 

0 = 

q^F3 

+ q^T3 = 

-g*Q 3 sin(r7 -b X2 -b Xg) - X3, 

(27) 

where |x||x| = to 

-2 and sgn(x); - x^)x(x|) = 

I are 


substituted according to the assumption: X 4 — Xg = 
X 4 = w > 0 and Xg > 0 in Section 3.2. 

Note that the equations (25), (26), and (27) represent 
the balance of the torque from the brake force and the 
gravity force at about 6 * 1 , 6 * 2 , and 63 , respectively. In 
particular, (25) can also be derived from the balance 
of the energy supply from the gravitational potential 
and the energy consumption via Coulomb friction and 
quadratic resistance. 

The steady-state equations in (24), (25), (26), and 
(27) can be reduced to the following form: 


r x; = X* = 0 , 

v .*\2 


ci{xX) = sin(a) —/x*X 3 > 0 , 


sin(x;) = - 
sm(x 2 + Xg 


9*Qi 
■f ry) = - 


9*Q2, 


(28) 


Thus, we have derived the steady-state equations of 
the FCWIP model with three unknowns X 2 , Xg, and x| 
(= uj). Note that the angle of the wheel in the steady 
state X* (t) oc ut never appears explicitly in these steady- 
state equations. 
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It is clear from (28) that under the given mechanical 
structure and environment, the nontrivial components 
of steady state X 3 , a;|) depend on three parameters, 

namely /i*, 77 , and ci. More precisely, and can 
be solved independently of x\, and they depend on the 
nondimensional friction /r* and the offset rj of the FCM. 
After that, x\ is obtained as a function of depending 
on Cl- 


3.4 Jacobian matrix at steady state 

We consider a variation 6 x of x around x* as x := 
X* + 5x and substitute it into the state-space model 
(16) as M.{x* + (5a;){a;* -I- <5®} = (/ -|- t)(£c* -I- 6 x). The 
zth component of the left side can be written by the 
Einstein notation as 

Li = Mij{x* + (5£c)|i:* -f 

= ^Mij{x*) -f -|- 0 {6x, Sxj)^ 

= Mij{x*)x* -\- Mtj(x*)Sxj 

+ ^^^Sxkx] + 0(Sx,Sx)^. ( 29 ) 

OX}^ 

Due to the structures of A4 = diag(E 3 ,M) and x* = 
{oj, 0,0,0,0,0)^, we have 


dM^j 

dxk 


SxkX* 


= 0 


0 • dxkUJ 


dM 


dxk 


—Sxk ■ 0 


(j = 1) 
(j > 1) 


(30) 


for all i and j. Therefore, neglecting the second and 
higher order term of Sx and 6x, we arrive at a variation 
equation of (16) as 


6x = M{x*) ^^Df{x*) -f Dt{x*)'^6x =: J{x*)Sx, 

(31) 


where := cos(x 2 ), C| := cos(a ;2 + x^ +ri), and 


Dt{x*) 



0 

0 

0 

0 

0 

\ 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 


—2cix| 0 - 

-Cbf^ 

N 

0 

0 

* 

0 

0 

* * 


Vo 

0 

-1 

0 

0 

— r* 
‘-b 

/ 


(33) 


To obtain (32) and (33), |a; 4 |a ;4 = x1 and sgn(x 4 — 
X 5 )x(x 3 ) = 1 are assumed, because xs(f) > 0 , X 4 {t) > 0 , 
and X 4 {t) — x^{t) > 0 hold for Xi{t) = x* -\- Sxi{t), 
|fei(f)| <1 (i = 1, • • • ,6). 

The results in (32) and (33) imply that the stability 
of the FCWIP model depends on the nondimensional 
viscous coefficient cl in addition to the parameters /x*, 
rj, and ci that determine the steady states (xl,,xl,xX). 

In summary, we have found that under a given me¬ 
chanical structure and environment, 


— the steady angles (xl^xD depend on {fj,*,r]), 

~ the steady descent velocity a;| depends on (/x*, ? 7 , ci), 
and 

— the stability depends on (/x*, 77 ,ci,c^). 


3.5 Eigenvalue equation 

It can be proved that rank [J(a;*)] = 5 < 6 = dim a:, 
which follows from the assumption of the uniform motion 
xi{t) = OJ. Thus, the characteristic equation of J{x*) is 
given in the following form: 

det(J(a;*) - AEe) = Ah(A) = 0, (34) 

/ 7 (A) = A^ -f UiA^ -j- U 2 A^ UsA^ -f (X 4 \ -\- (X 5 , (1^^) 

where det(-) denotes a determinant and Eg is a 6 x 6 
identity matrix. Eor simplicity, we refer to h(X) = 0 as 
an eigenvalue equation of the FCWIP model. 

Therefore, the steady state x* becomes stable if the 
maximal real part of the eigenvalues is negative, that is 


where Df{x*) denotes the Jacobian matrix of f{x) 
around x* and J(a:*) provides a closed-loop state matrix 
whose eigenvalues represent the stabilities of the FCWIP 
model. The components of J{x*) are given by 

0 0 I 0 0\ 

0 0 0 1 0 

0 0 0 0 1 

0 0 000 ’ 

I 0 g*Q4Cl - g*QsCi -g*Q^Cl 0 0 0 
Vo -g*QiCl -g*QzCl 0 0 0/ 


Df{x*) = 


/o 

0 

0 

0 


(32) 


A := Re Amax < 0 , A^ax := argmaxRe (A), (36) 

A 

where Ai (z = 1, • • • ,5) are the roots of h(A) = 0. 


4 Critical points of steady state 

As mentioned in Section 3.4, the stabilities of the steady 
state X* depend on /x*, 77 , ci, and cj. Here, we provide 
some numerical examples of that dependency. Thus, the 
parameter values are set to those listed in Table 1, 2, 
and 3 by default unless otherwise noted. 
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Fig. 5 Maximal real part of eigenvalue A and steady state 
X* as functions of the nondimensional friction /j.* . 



4.1 Dependency on /i* and y 


We first examine the dependency on the nondimensional 
friction fx* and the offset r/, which determine the steady 
angles and Xg via (28). 

Figure 5 shows the maximal real part of the eigen¬ 
value A and the nontrivial components x* {i = 2, 3, 4) of 
the steady state x* as functions of the nondimensional 
friction fi* . A and x* are obtained from numerical solu¬ 
tions of (28) by Newton’s method and (36), respectively. 
The solid line in the top graph represents A, and the 
solid and dotted lines in the bottom graph represent 
stable and unstable x*, respectively. The values of x* 
{i = 2 ,3) are scaled to share a common vertical axis. 

It is clear from Fig. 5 that three types of critical 
points Pq, Pi, and P 2 appear, which are denoted as 
open circles, filled circles, and a triangle, respectively. 
Pq gives an infimum inf^* (x))) = 0 of the descent velocity 
3^4 > 0, Pi gives a stability boundary, and P 2 gives a 
folded (nonsmooth) minimum of A{p,*). 

To characterize these points. Fig. 6 plots the root 
loci of h{X) = 0 in (35) along the steady states x* 
in Fig. 5. It can be numerically proven that under the 
considered condition, the eigenvalue equation h{X) = 0 
has three real roots and a pair of complex conjugate 




X 3 (xlO^) 

Fig. 7 Limit cycle for p* = 1.01 (just after Hopf bifurcation). 


roots as 

Aq Sq, Ai 4 ^ Si ib jI2, A 2 S 2 , Ag Sg, 

So, Si > S 2 ,S 3 < 0, (37) 

where Sj, 17 are real numbers and j := 3 /—!. Among the 
five roots, only the first three roots Aq, Ai+, and Ai_ 
affect the stability change, and only these three roots 
appear in the range of Fig. 6 . 

On the basis of the root loci obtained in Fig. 6 , we 
can characterize the critical points in terms of eigenvalue 
types as follows: 

— (PO): Pq is a point such that Amax = 0. 

— (PI): Pi is a Hopf bifurcation point, for which the 
loci cross the imaginary axis at Amax = ijl7, where 
17 is an angular frequency of a limit cycle. 

— (P2): P 2 is a point such that maximal real root sq 
and the real part of the complex conjugate roots Si 
coincide, at which point they switch roles to produce 
the maximal real part. 

Physically speaking, Pq and Pi provide stability limits, 
and P 2 maximizes the total stability or minimizes the 
time constant of the FCWIP model. The descent velocity 
vanishes (X 4 = 0) at Pq in this case, and a self-excited 
oscillation (limit cycle) emerges for A{fx*) > 0 near Pi. 

Figure 7 (a) shows the limit cycle for /i* = 1.01 
immediately after the Hopf bifurcation point Pi. Under 
this condition, the FCWIP model descends the slope in 
a standing position while the angles of the pendulums 
oscillate slightly and periodically. As this limit cycle is 
stable. Pi is identified as a Hopf bifurcation point. In 

























Katsutoshi Yoshida et al. 


general, a limit cycle occurs because of the temporal 
presence of negative resistance (i.e., negative energy 
consumption per unit time). In our problem, this is 
given by the braking torque T 2 in ( 20 ) multiplied by 
the friction velocity (xa — X 5 ), specifically as 

V = T 2 ■ (X 4 - X 5 ) 

= q~^fj,*sgn(x4 - X5)x(x3)(x3 + clxe) ■ {X4 - X 5 ) 

= q~‘^fj^*\x4 - X5\x{x3){x3 + clxe) 

=: C'x(a; 3 )(a ;3 + cIxq) (C > 0), (38) 

where sgn(x)a; = \x\ is applied. As x( 2 ^ 3 ) is a step 
function, 22 < 0 occurs when a ;3 > 0 and X 3 + c^xq < 0 . 
This implies that when 22 < 0, the brake pad (see Fig. 2 
recalling z = px 3 ,z = {p/q)xe) will be released with a 
velocity less than xq < —X 3 /CI < 0 while the pad is still 
pressed against the disk {xs > 0). Figure 7 (b) shows 
the energy consumption 22 along the limit cycle as a 
function of X 3 , which numerically confirms the presence 
of 22 < 0. During this negative energy consumption, the 
cycle comes across the deadband border 3:3 = z = 0. At 
this point, the step function x(a: 3 ) jumps from 1 to 0 . 
This causes the 22 to jump from a negative to a positive 
value, similar to xq. 

Figure 8 (a) shows the result as functions of the offset 
rj for p* = 0.97. In this case, a Hopf bifurcation point 
Pi does not appear in the plotted range 77 > 0. Outside 
this range, X 2 and x^ vanish at 77 = 0 and violate the 
physical assumptions < 0 and > 0 for 77 < 0 . 


4.2 Dependency on ci and cj 

Figure 8 (b) shows the results as functions of the quadratic 
resistance ci. It is clear that all the critical points Pq, 
Pi, and P 2 appear, although X 2 and X 3 become constant 
here because ci only affects x\, as already discussed in 
(28). 

However, the physical results of Pq are different. That 
is, Pq (or Amax = 0) On A{p*) in Fig. 5 corresponds to 
the descent velocity at rest = 0. In contrast, Pq on 
A(ci) in Fig. 8 (b) corresponds to the infinite descent 
velocity a;| —>■ oo (ci —>■ 0). This can be explained by 
the second equation in (28), which is hyperbolic with 
respect to ci and x\'. 

ci{x\Y = g*Qism{a) — p*x"^ =:C>Q. (39) 

This equation exhibits the following features: 




0 2 4 6 8 10 

Cj {xlO“^) 



1 23456789 10 


Fig. 8 Maximal real part of eigenvalue A and steady states 
X* for p* = 0.97 as functions of (a) the offset 77 , (b) the 
quadratic resistance ci, and (c) the viscous coefficient . 


The second feature ci(xl)'^ = 0 holds when ci = 0 
and/or xl = 0. The latter condition 0:4 = 0 directly 
explains Pq in Fig. 5 for a constant ci > 0. Onthe other 
hand, Pq in Fig. 8 (b) can be explained by the limit 
Cl — >■ 0 that causes xl = \JCjci -A 00 for the constant 
C. In addition, these different Pq can also be explained 
physically. That is, the condition 01 ( 0 : 4 )^ = 0 results in 
vanishing of the quadratic resistance force ci\xWx\ = 0 . 
This can be caused by the mechanism at rest x\ = Q 
as well as by the absence of the effect of the quadratic 
resistance ci = 0 . 


— The right side of (39) (e.g., C) is expected to be 
constant because Xq is determined independently of 
(39). 

— The left side ci{xl)‘^ vanishes at Pq (or Amax = 0), 
as will be discussed in Section 5.1.1. 


Figure 8 (c) plots the result for the nondimensional 
viscous coefficient cl of the FCM. It is found that only the 
Hopf bifurcation point Pi appears. Moreover, it appears 
that X* {i = 1,2,3) are all constant because the steady- 
state equation in (28) is independent of cl, which only 
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affects the components of the Jacobian matrix as 
and —cl in (33). 


5 Stability limits 

Finally, we numerically continue the critical points in 
two-parameter planes to characterize the stability limits 
of the FCWIP model. 


5.1 Conditions of the critical points 

5.1.1 Zero quadratic resistance 

The condition of Pq (or Amax = 0) is mathematically 
equivalent to 05 = 0 in the eigenvalue equation (35). It 
follows that 


eliminating 17 from them, we obtain 

:= a5|a4(—010203 -h O3 -h 0^04) 

“h ( — 02 O 3 -|- 01(02 — 204 )) 05 -|- Ogj’ =0 (42) 

for the Hopf bifurcation point. It is implied from (42) 
that Pi coincides with Pg (or P§°) because 05 = 0 for 
Pq (or Pg°°) leads to <?i = 0. 

Note that, rigorously speaking, the above condition 
provides only a necessary condition of a Hopf bifurcation 
point; however, it leads to satisfactory results in the 
present analysis. 

5.1.3 Minimal time constant 

We numerically detect the condition of P 2 for the min¬ 
imal time constant that satisfies 


0 = 05 = AeiM ^ ■ 2cix\g* 

X { - ClQi - C*MClg*Qi - (1 + /X*)}} 
^ 0 = cixl[ciQi + C^Qii{Clg*Qi - (1 + /x*)}} 
cia ;4 = 0 

Cl = 0 or a ;4 = 0 ci{xXY =Q. (40) 

Therefore, it is mathematically shown that the sufficient 
condition for Pg (or Amax = 0) is given by the zero 
quadratic resistance 01 (^ 4 )^ = 0 . 

As discussed in Section 4.2, the condition Pg (or 
Cl ( 2 : 4 )^ = 0 ) causes two distinct descent velocities: xX = 
0 for Cl > 0 and = 00 for ci = 0. Therefore, we 
denote them as Pg and P“, respectively, in the following 
sections. 

The condition of Pg is given by (39) with ci ( 2 : 4 )^ = 0, 
from which we can eliminate x^ and Xg to obtain 


|so-si|<10-® (43) 

where sg is the maximal real eigenvalue and si is the real 
part of the complex conjugate eigenvalues, as defined 
in (37). In numerical calculations, the parameter values 
considered are swept to detect a point that satisfies 
(43), where the point in the first detection is taken as 
the point detected. 

Note that we attempted to derive an equation for 
the minimal time constant in a closed form of oi, • • • ,05 
based on a given form of eigenvalue equation: h{X) = 
(A Amax)(A Amax Jc)(A Amax T 4“ pA J- q) ^ 

however, the result was very weak to detect P 2 pre¬ 
cisely. Therefore, in this paper, we employ the numerical 
method mentioned above, although another approach 
would be possible for an analytical expression of P 2 . 

5.2 Numerical continuation of the critical points 


:= X 7 + arcsin 


QisinaA g*(5isina 


— arcsin 


Q 3 / 2 * / /X 

Qi{l + fJ,*) sin a 


Qifi* 


= 0. (41) 


On the other hand, P“ is simply given by ci = 0. 
Because the condition (41) does not contain x^ and 2 : 3 , 
the zero quadratic resistance points Pg and P§° are 
determined independently of the pendulum angles x^ 
and x^. 


5.1.2 Hopf bifurcation 

The condition Pi (or Amax = ±jI2) for a Hopf bifur¬ 
cation point is given by Keh{jf2) = lmh{jf2) = 0. On 


Figure 9 plots the sets of the critical points Pg , Pi, and 
P 2 on two-parameter planes obtained from the numerical 
solutions of (41), (42), and (43) under (28). The results 
on the (/x*,? 7 ), (/x*,ci), and (/x*,Cj) planes are labeled 
(a), (b), and (c) respectively in Fig. 9. The solid line 
denotes the set of the zero quadratic resistance point Pg 
for X 4 = 0, the dotted line denotes the set of the Hopf 
bifurcation point Pi, and the chained line denotes the 
set of the minimal time constant point P 2 . The plots 
of P§° for Cl = 0 do not appear in the ranges of these 
plots. The hatched areas represent the stable regions of 
the steady state satisfying A < 0 in (36). 

It is clearly seen in Fig. 9 that the stable regions 
are bounded by Pg and Pi and that the minimal time 
constant point P 2 is sandwiched between them. Note 




10 


Katsutoshi Yoshida et al. 





n* 

Fig. 9 Continuation of the critical points Pq, Pi, and P 2 
on the two-parameter planes: (a) for the (/r*, rj) plane, (b) 
for the (/r*,ci) plane, and (c) for the plane. The 

hatched areas represent asymptotically stable conditions of 
the steady states. 


that in Fig. 9 (a), the plots are bounded by the assump¬ 
tion rj > 0 and that in Fig. 9 (b) and (c), Pq lies along 
the vertical line at ^* = fi* ~ 0.89474. This is because 
Pq is determined independently of ci and in (41). 
Moreover, it appears that the critical points Pq, Pi, and 
P 2 tend to coincide as fj,* increases on the (/i*, 77 ) plane, 
as cl increases and decreases on the {^*, cl) plane, and 
as Cl increases on the (/t*,ci) plane. In contrast, as ci 
decreases on the (/i*,ci) plane, Pi and P 2 also tend 
to coincide but they approach ci = 0 or P“ instead 
of Pq. As discussed in Section 5.1.2, the convergence 
between Pq (or P“) and Pi can be explained by (42), 
in which the condition ^1 = 0 for the Hopf bifurcation 
point Pi contains 05 = 0 for the zero quadratic resis¬ 
tance point Pq (or P§°). As shown in Fig. 10, however, 
it can be numerically proven that a purely imaginary 
eigenvalue Amax = along Pi does not vanish even 
when Pi coincides with Pq (or P^) in the parameter 
planes. Therefore, the Hopf bifurcation point Pi does 


0.45 
0.4 
0.35 
^ 0.3 

0.25 

0.9 1 1.1 1.2 1.3 1.4 

fi* 

Fig. 10 Imaginary part of complex conjugate eigenvalues 
along Pi in Fig. 9 (a). 



not degenerate and double-zero eigenvalues never arise 
at that point. 

From an engineering viewpoint, the results obtained 
above imply that the FCWIP model allows a certain 
amount of tolerance in parameter settings because the 
stable conditions are obtained as simply connected finite 
areas. This suggests the possibility that the proposed 
mechanism can work even if there are some manufactur¬ 
ing errors. Furthermore, the stable area on the {fj,*,r]) 
plane in Fig. 9 (a) forms a monotonically increasing 
shape, which implies that the offset rj can be designed 
to shift the stable range of the friction coefficient n*. 


5.3 Numerical evaluation of the descent velocity 

In view of engineering applications, the descent velocity 
(or angular velocity of the wheel) a;| must be adjusted 
to a value suitable for the intended use. Figure 11 shows 
the values of a;| mapped into a gray scale within the 
stable areas on the parameter planes in Fig. 9. values 
are numerically obtained by solving (28). 

It is clear from Fig. 11 that tends toward zero 
as the parameter conditions approach Pq (solid lines), 
which is in agreement with the definition of Pq. Es¬ 
pecially, in Fig. 11 (b), it is also clarified that a;| is 
diverging as the condition approaches P“ or ci = 0 . 
Moreover, it appears that changes smoothly and 
monotonically with the variations of the parameters. 
This suggests that one can adjust the descent velocity 

by continuously shifting the parameters. However, it 
is also observed in Fig. 11 (b) that a sufficiently large 
value of the quadratic resistance ci is required to sta¬ 
bilize the mechanism because although a range of /i* 
exists for small a;| near /i* = /i* « 0.89474, it narrows 
significantly as ci decreases. In addition, it appears in 
Fig. II (c) that X 4 does not depend on cl, as discussed 
in the last paragraph in Section 3.4. 
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Fig. 11 Angular velocities of the wheel 2:4 within the stable 
areas on the two-parameter planes: (a) for the (/r*,?;) plane; 
(b) for the (/i*,ci) plane; and (c) for the plane. 


6 Conclusion 

For a non-electrified method to stabilize a wheeled in¬ 
verted pendulum descending a slope, we proposed the 
FCWIP mechanism, a wheeled double pendulum, whose 
second pendulum transforms gravity force into brake 
force acting on the wheel. We conducted steady-state 
analysis of the proposed model and obtained the fol¬ 
lowing results: 

— The steady angles depend on 

— The steady descent velocity depends on (/i*, 77 , ci). 

— The stability depends on (/i*, 77 ,ci,c^). 

Then, we found three types of critical points in the 
steady states, as 

— Pq: the point for zero quadratic resistance on the 
wheel. (Pq for X 4 = 0 and for X 4 = 00 in detail) 

— Pi: the Hopf bifurcation point. 

— P 2 : the point for the minimal time constant. 


Finally, we conducted numerical continuations of 
these points on the two-parameter planes and evaluated 
the descent velocity to obtain the following results: 

— The stable conditions are obtained as simply con¬ 
nected finite areas on the parameter planes, bounded 
by Pq and Pi. 

— The minimal time constant point P 2 is sandwiched 
between Pq and Pi. 

— The descent velocity x\ changes smoothly and mono- 
tonically with the parameter variations. 

The abovementioned results lead to the conclusion 
that the parameter selection to design the FCWIP mech¬ 
anism stabilized on a slope will not be highly sensitive, 
at least in theory. 

In future work, we plan to develop a physical FCWIP 
mechanism. For this purpose, we will introduce stick- 
slip effects into the friction term in our model and 
investigate the effect on the stabilities. We also plan 
to conduct stochastic analysis on the FCWIP model to 
consider robustness against random disturbances and 
random parameter fluctuations. 
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